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ABSTRACT 

Ionizing UV radiation and supernova flows amidst clustered minihalos at high redshift regulated 
the rise of the first stellar populations in the universe. Previous studies have addressed the effects 
of very massive primordial stars on the collapse of nearby halos into new stars, but the absence 
of the odd-even nucleosynthetic signature of pair-instability supernovae in ancient metal-poor stars 
suggests that Population III stars may have been less than 100 Mq. We extend our earlier survey of 
local UV feedback on star formation to 25 - 80 M stars and include kinetic feedback by supernovae 
for 25 - 40 M Q stars. We find radiative feedback to be relatively uniform over this mass range, 
primarily because the larger fluxes of more massive stars are offset by their shorter lifetimes. Our 
models demonstrate that prior to the rise of global UV backgrounds, Lyman- Werner photons from 
nearby stars cannot prevent halos from forming new stars. These calculations also reveal that violent 
dynamical instabilities can erupt in the UV radiation front enveloping a primordial halo but that 
they ultimately have no effect on the formation of a star. Finally, our simulations suggest that relic 
H II regions surrounding partially evaporated halos may expel Lyman- Werner backgrounds at lower 
rcdshifts, allowing stars to form that were previously suppressed. We provide fits to radiative and 
kinetic feedback on star formation for use in both semianalytic models and numerical simulations. 
Subject headings: cosmology: theory — early universe — galaxies: high redshift — H II regions — 
intcrgalactic medium — radiative transfer 



1. INTRODUCTION 

The survival of cosmological minihalos in the ionizing 
and Lyman- Werner (LW) UV fields of primordial stars 
is key to the rise of stellar populations at high redshifts. 
Numerical models suggest that the first stars are very 
massive, 25 - 500 M Q , and that they form in isolation in 
small dark matter halos of ~ 1 5 - 10 7 M w at z ~ 20 
- 30 (IBromm et all 119991 120021: INakamura fe Umemural 
l200lt lAbel et all 120001 120021: IQ'Shea fc Normanl 12001 7 
These stars create large H II regio ns 2.5-5 kpc in ra- 
dius that can engulf nearby halos (IWhalen et all [20041: 
Kitayama et all 120041: lAlyarez et all 120061: lAbel et all 
20071 IWise fc Abell 12008^ From z ~ 10 - 20 Pop- 
ulation III (Pop III) stars also build up a global LW 
background that sterilizes minihalos of H2, delaying or 



Ahn fc Shapirdl2"0"07t IWhalen et al|[2008al: IWise fc Abell 
20081* iHaseeawa et al.H2009t ISusa et al.ll2009D and with- 
out (|Shapiro et al.ll2004t llliev et al.ll2005D H? gas chem- 
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20091 lAhn et alJl2009t l . At high redshifts. ionizing radi- 
ation is relatively local while LW photons can originate 
from many megaparsecs away because their energies lie 
below the ionization limit of H. 

The picture is simpler in the first generation, in 
which there is no LW background. Simulations 
of halo photoevaporation by nearby very massive 
(> 100 M (?)) primordial star s have been performed , 
both with (IQ'Shea et all [20051 : ISusa fc Umemural [20061 : 
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istry. The two- and three-dimensional studies are in 
good agreement with each other but not w ith the one- 
dimensional studies (|Ahn fc Shapiro! 120071 ). primarily 
due to the serious hydrodynamical artifacts that arise 
from unphysical shock focusing and central bounce that 
occur in one-dimensional spherical coordinate meshes. 
When multifrequency photon transport is coupled to pri- 
mordial chemistry and hydrodynamics, halo photoevap- 
oration proceeds in two stages. First, LW photons from 
the star reach the halo before the I-front, partly disso- 
ciating molecular hydrogen in its core and slowing or 
halting its collapse. Not all of the H 2 is erradicated be- 
cause molecular hydrogen deep in the halo partly shields 
itself from the photons and small free electron fractions 
restore H2 by catalysis through the H~ channel. The 
I-front decelerates as it approaches the halo and trans- 
forms from R-type to D-type. Molecular hydrogen forms 
in the small ionized fractions and warm temperatures 
in the outer layers of the front (e.g. iRicotti et al.l 120011 : 
IWhalen fc Normanll2008al ). partly shielding the interior 
of the halo from LW photons from the star and even 
allowing H2 fractions in the core to rise in some cases. 
The ionization front preferentially advances through the 
low densities above and below the halo and assumes a 
cometary shape, with the inner regions of the halo cast- 
ing a shadow. 

Second, after the death of the star the H II re- 
gion surrounding the halo begins to recombine out 
of equilibrium, rapidly forming H2 and HD that can 
cool the ionized gas down to the CMB temperature 
(INakamura fc Umemural 120021: Johnson fc Bromml 120061 : 
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H2 reconstitutes in the core, rapidly surpassing its orig- 
inal levels and allowing gas to again cool and collapse. 
Meanwhile, the I-front shock remnant continues to con- 
verge on the halo core from one side while warm ionized 
gas presses down into the shadow and wraps around the 
halo from behind. If the UV flux from the star is not too 
high, the relic shock compresses the core and enriches it 
with the H2 originally formed in the I-front, accelerat- 
ing its cooling and collapse into a new star. If the flux 
from the star is large, the ionized gas instead disrupts the 
core and pre-empts new star formation. The studies per- 
formed to date assume stars that do not die in supernova 
(SN) explosions and therefore exc lude ram pressure strip- 
ping of the halo by th e remnant (|Cen fc R iquclme1 l2008t 
iSakuma &: Susall2009l ) and its contamination and cooling 
by heavy elements. In the first generation of stars, both 
stages of halo photoevaporation occur on time scales that 
are short in comparison to merger or accretion times. 

In large scale calculations of cosmological structure for- 
mation, halo evaporation is usually modeled with meta- 
galactic ionizing and LW backgrounds that uniformly 
permeate the simulation volume. Such halos evolve very 
differently than when photon transport is performed be- 
cause they are photoevaporated and photodissociated 
from the inside out. Rather than being compressed and 
shielded from LW photons, baryons are expelled from all 
directions. Consequently, radiative feedback is invariably 
negative in these mod els. An excepti o n to t his are the 
recent simulations by iMcsingcr et al.1 (p009^ . who find 
that low mass halos that are easily evaporated at early 
times later pool gas and form stars when the uniform 
fields are turned off and the fossil H II regions cool down. 
Self-shielding of H2 from external UV sources cannot be 
modeled by imposed backgrounds, so molecular hydro- 
gen is driven to much lower levels in these halos than in 
real ones. In general, negative feedback in studies that 
rely on metagalactic UV backgrounds should be taken 
as a (possibly extreme) upper limit. Star formation in 
the UV environments of the early universe is likely to be 
significantly more robust. 

In principle, the parameter space for local radiative 
feedback between clustered minihalos is vast even if su- 
pernova effects are excluded, making the formulation of 
general rules for use in analytical models problematic. 
Halo mass, central baryon density, distance and lumi- 
nosity of the star, stellar lifetime and spectral profile all 
govern halo evaporation. Distance and luminosity can be 
combined in a single parameter over intervals in stellar 
mass for whi ch the shape of the sou rce spectrum does not 
vary much. iWhalen et al.l (|2008al) further reduced this 
parameter space by examining feedback in the smallest 
halo ever found to form a star in an adaptive mesh refine- 
ment (AMR) calculation but at four stages of collapse. 
Since more massive halos at the same central baryon den- 
sity would be less affected by radiation, their findings 
constitute a conservative upper limit to the influence one 
star can exert on any neighbor halo capable of forming 
a star. 

In this paper we extend our previous survey by exam- 
ining the evaporation of minihalos by low-mass primor- 
dial stars prior to the buildup of a large LW background. 
Pop III stars from 25 - 80 M Q have spectra with larger 
LW/ionizing UV ratios than those of the 120 M Q star in 
our previous study and they illuminate other halos for 



longer times. We sample the same halo and central gas 
densities as in our earlier study in order to place upper 
limits on the radiative feedback of such stars on more 
massive halos. By determining the final state of the halo 
core several Myr after the death of the star we construct 
rules for local radiative and kinetic feedback as a function 
of stellar mass, initial core density, and distance to the 
star. In § 2 we review our numerical methods, in § 3 we 
describe our cosmological halo models, in § 4 we tabulate 
both radiative and kinetic feedback in a photoevaporated 
halo, and in § 5 we conclude. 

2. NUMERICAL ALGORITHM 

We perform o ur halo photoeva poration calculations 
with ZEUS-MP (jHaves et al.ll2006f) . a massively-parallel 
Eulerian astrophysical hydrodynamics code that solves 
the equations of ideal fluid dynamic^: 
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Here, p, e, and the Vi are the mass density, inter- 
nal energy density, and velocity at each mesh point 
and p — (7 — l)e and Q are the gas pressure and 
the von Neumann-Rich tmeyer artificial viscosity tensor 
(jStone fc NormanlH992T) . ZEUS-MP evolves these equa- 
tions wit h a second-orde r accurate monotonic advection 
scheme (jvan Leerl 119771 ) in one, two, or three dimen- 
sions on Cartesian (XYZ), cylindrical (ZRP), or spheri- 
cal polar (RTP) coordinate meshes. Our augmented ver- 
sion of the publically- available code self-consistently cou- 
ples pr imordial gas chemistry ()Whalen fc Normanil2006l 
2008b) and multifrequency photon-cons erving UV radia- 
tive transfer (Whalen fc Normanll2008al ) to fluid dynam- 
ics for radiation hydroynamical transport of cosmological 
I-fronts. 



2.1. Primordial H and He Chemistry 

We evolve H, H+, He, He+, Hc 2 +, H 
e~ with nine additional contin uity equations and t he 
nonequilibrium rate equations of lAnninos et al.l (|1997f) : 



Hj , H 2 , and 



^ = - V ■ ( Pi w) Pih(X)PiPh - E Wii ( 4 ) 



where (3jk is the rate coefficient of the reaction between 
species j and k that creates (+) or destroys (-) species 
i, and the Kj are the radiative rate coefficients. We as- 
sume that the species share a common velocity distri- 
bution. Mass and charge conservation, which are not 
guaranteed by either chemical or advective updates, are 
enforced each time the fluid equations are solved. The 
divergence terms and reaction network are operator-split 
and evolved on their respective time scales, as we explain 
in greater detail below. 

5 http:/ /lea. ucsd.edu/portal/codes/zeusmp2 
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Microphysical heating and cooling due to photoioniza- 
tion and gas chemistry is coupled to the gas energy den- 
sity by an isochoric update that is operator-split from 
updates to the fluid equations: 



= r-A, 



(5) 



where T is the cumulative heating rate due to photons 
of all frequencies and A is the sum of the cooling rates 
due to collisional ionization and excitation of H and He, 
recombinations of H and He, inverse Compton scatter- 
ing (IC) off the CMB, brem sstrahlung emission, and H2 
cooling (|Galli fc Pallalll998h . 

2.2. Radiative Transfer 

Ou r photon-conserving UV transport (|Abel et al.l 
119991 : iMellema et~aH 120061 ) . which is distinct from the 
flux-limited diffusion native to the public release of 
ZEUS-MP, solves the static approximation to the equa- 
tion of transfer in flux form to compute radiative rate 
coefficients for the r eaction network at every point on 
the coordinate mesh (|Abel et al.lll999f) . As currently im- 
plemented, our code can transport photons from a point 
source centered in a spherical grid or in plane waves along 
the x or z-axes of Cartesian or cylindrical meshes. In our 
models, the radiation from the star is treated as a plane 
wave that propagates along the positive z-axis. Plane 
waves are a good approximation to the UV ionizing flux 
at the interhalo separations in our study, but we atten- 
uate their intensity by 1/R 2 to account for geometrical 
dilution. 

As in iWhalen et~aT1 (|2008aD . we discretizc the black- 
body photon emission rates of the stars in our survey 
with 40 uniform bins from 0.755 to 13.6 eV and 80 log- 
arithmically spaced bins from 13.6 eV to 90 eV, again 
normalizing the m by the total ionizing photon rates for 
Pop III stars bv lSchaererl (|2002j ). Th e radiative reactions 
in our m odels are listed in Table 1 of lWhalen fc Normanl 
( 2008al). We do not evaluate H 2 photodissociation rates 
with radiative transfer. Instead, we calculate them along 
rays parallel to the direction of radiation flow using self- 
shielding func tions modified for t herm al broadening as 
prescribed by iDraine fc Bcrtoldi (|1996l ) to approximate 
the effects of gas motion. They are shown in equations 
9 and 10 of lWhalen fc Normanl (|2008aD . 

2.3. Radiation Forces 

Since our prior survey of radiative feedback, we have 
implemented momentum deposition in the gas due to 
ionizations. Radiation pressure in ionizing UV transport 
comes into play at two locations: at the I-front itself and 
in reco mbining gas in the H II region. Wh alen fc Normanl 
(2006) examined the acceleration of fluid elements at the 
front and found that it was large but momentary, and 
that its inclusion alters the velocity of the front by only 
1-2 km s . Direct momentum deposition within the H 
II region is only prominent where gas is very dense, like 
the center of a cosmological minihalo being evaporated 
by a star at its center. There, rapid successive cycles of 
ionization and recombination can impart radiation forces 
to the gas that are hundreds of times the strength of 
gravity at early time s (lower left panel of Figure 1 in 
iKitavama et al.ll2004l ). As these forces propel gas near 
the center of the halo out into the H II region, its densities 



and recombination rates fall, so more ionizing photons 
from the star reach the I-front. This higher flux results 
in I-fronts that are faster than when such forces are not 
included. Thus, in early UV breakout, radiation forces 
should speed up the I-front, and in trial runs we find that 
D-type fronts are 10 - 20% faster than when momentum 
transfer due ionizations is neglected. However, this effect 
is transient: after the internal rearrangement of gas deep 
within the H II region dilutes its interior, radiation forces 
there sharply fall. 

We expect much weaker forces in this study because 
the front climbs a density gradient as it approaches the 
halo rather than descending one, so the gas behind the 
I-front is always relatively diffuse. Since new ionizations 
due to recombinations are far less frequent, less momen- 
tum will be imparted to the gas, so we include these 
effects only for completeness. Updates to the gas veloci- 
ties are straightforward since momentum deposition due 
to direct photons from the source is always parallel to the 
direction of radiation flow. We describe the time scales 
on which momentum updates are performed below. 

2.4. Adaptive Subcycling 

A hierarchy of highly disparate characteristic time 
scales arises when gas dynamics, radiative transfer, and 
primordial chemistry are solved in a given application. 
The three governing times are the Courant time, the 
chemical time 
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and the photoheating/cooling time 
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Their relative magnitudes can seamlessly evolve through- 
out a single application. For example, when an I-front 
propagates through a medium, photoheating times are 
often smaller than Courant times, and chemical time 
scales are usually shorter than either one. On the other 
hand, fossil H II regions can cool faster than they re- 
combine, so cooling times become shorter than chemi- 
cal times. The key to solving all three processes self- 
consistently is to evolve each on its own timescale with- 
out restricting the entire algorithm to the shortest of the 
times. To successfully deal with both I-fronts and relic 
H II regions, an algorithm must adaptive ly reshuffle the 
time scales on which the three processes are solved. Im- 
plicit schemes are sometimes applied to stiff sets of dif- 
ferential equations like those in our model because they 
are unconditionally stable over the Courant time. How- 
ever, accurate I-front transport in stratified media often 
requires restricting updates to both the gas energy and 
fluid equations to photoheating times in order to cap- 
ture the correct energy deposition into the gas, and lin- 
ear systems solves over such short time scales would be 
prohibitive in more than one or two dimensions. Enforc- 
ing photon conservation in implicit schemes can also be 
problematic. 

We instead subcycle chemical species and gas energy 
updates over the minimum of the chemical and heat- 
ing/cooling times until the larger of the two has been 
crossed, at which point we perform full hydrodynamical 
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Fig. 1. — Spherically-averaged baryon profiles for the f.35 X fO 5 Mq halo at four stages of evolution. The redshifts of the 023, 039, 059, 
and 073 profiles are 23.9, f7.7, 15.6, and 15.0, respectively, with corresponding central densities of 1.43, 10.5, 108, and 1596 cm -3 . Left: 
densities. Right: radial velocities. 



updates of gas densities, energies, and velocities. These 
times are global minima for the entire grid. The chemi- 
cal times arc dchncd in terms of electron flow to accom- 
modate all chemical processes rather than just ioniza- 
tions or recombinations. Adopting the minimum of the 
two times for chemistry and gas energy updates enforces 
accuracy in the reaction network when t c hem becomes 
greater than the, such as in relic H II regions. Our adap- 
tive subcycling scheme is d escribed in greater detail in 
(|Whalen fc No rma n1 l2008blh 

3. HALO MODELS 

The 1.35 x 10 5 M Q halos we study are computed 
from cosmological initial conditions in the Enzo adap- 
tive m esh refinement (AMR ) code as described in de- 
tail in iWhalen et al.l ()2008aD . The halo baryon density 
and temperature profiles were spherically averaged and 
then mapped onto a two-dimensional cylindrical coordi- 
nate (ZR) grid in ZEUS-MP. The assumption of radial 
symmetry in two dimensions is approximate but valid, 
given the spheroidal nature of the halos. A single, sym- 
metrized baryon profile better represents all halos of this 
mass than the three-dimensional profile of just this halo 
because small differences in angular distribution of the 
baryons are averaged out. We center the halo on the 
z— axis so that only its upper hemisphere resides on the 
grid. The mesh boundaries are -125 pc and 125 pc in 
z and 0.01 pc and 125 pc in r. The grid is 1000 zones 
in z and 500 zones in r with a spatial resolution of 0.25 
pc. Outflow conditions are assigned to the upper and 
lower boundaries in z and reflecting and outflow condi- 
tions are imposed on the inner and outer boundaries in 
r, respectively. 

Dark matter gravity is included by constructing a sep- 
arate potential that cancels pressure forces everywhere 



on the mesh, thus imposing hydrostatic equilibrium on 
the halo, and then holding the potential fixed through- 
out the simulation. Force updates to gas velocities are 
performed with the potential every hydrodynamical time 
step during the simulation. Neglecting dark matter dy- 
namics introduces no serious errors because the gas in 
the halo evolves on much shorter time scales than cither 
the Hubble time or merger times, which at z ~ 20 are 
approximately 20 Myr. Updates to the self-gravity of the 
gas are also performed every hydrodynamical time step 
by evolving Poisson's equation with a two-dimensional 
conjugate gradient (CG) solver. 

In reality, the halo does have infall velocities as shown 
in Fig[TJ but they are minor. Their inclusion would only 
slightly enhance core densities during photoevaporation. 
The composition of the gas in our models was primordial, 
76% H and 24% He by mass. We assume ionized and H 2 
fractions of 1.0 x 10~ 4 and 2 x 10~ 6 , respectively. These 
values are consistent with remnant free electron fractions 
from the era of recombination at the redshift z = 20 
we take for our models and with those expected from 
collisional ionizations in cosmological accretion shocks. 
We consider the four evolutionary stages of the halo in 
our earlier work, corresponding to central densities n c = 
1.43, 10.5, 108, and 1596 cm, , respectively. Each halo 
profile is illuminated at 150, 250, 500, and 1000 pc, which 
are typical of inter halo separations in a cluster. Each 
profile is illuminated by 25, 40, 60 and 80 M Q stars for 
their respective lifetimes, 6.46, 3.86, 3.46 and 3.01 Myr, 
for a total of 64 models. These stars sample the lower end 
of the Population III mass spectrum. After the death of 
the star the halo is further evolved to a total of 10 Myr, 
or about half of a merger time, to determine the effect of 
relic H II region dynamics on the gas in the core. 
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3.1. Criteria for Star Formation 

Unlike AMR or Lagrangian SPH algorithms, the fixed 
Eulcrian grid in ZEUS-MP lacks the resolution to fol- 
low the collapse of the baryons in the halo into a new 10 : 
star. However, numerical simulations demonstrate that 
primordial star formation is a very robust process: if a 
halo is in a state in which its baryons can form a star, 
it will within a few Myr. Since our halo is one in which 
star formation has already been proven to occur, if the 
final density and H2 fraction of its core is greater than 
or equal to its pre-illumination levels, a star will form at 
its center. In some cases the relic I-front shock passes 
through the center of the halo and imparts an impulse to 
the gas there. However, if the final velocity of the per- 
turbed baryons within a parsec of the core is less than the 
escape speed from the halo, they will be trapped by the 
dark matter potential and still form a star. We therefore 
adopt the simple criteria that if central baryon densities 
and H2 fractions in the halo are the same or greater than 
those just before illumination, and if central gas veloci- 
ties are lower than the escape speed, a star will form in 1 0' 
the halo. If we approximate the binding energy E% of the 
baryons to the halo by that of a homogeneous sphere, 
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where Mh, nib, and i? v ir are the total halo mass, baryon 
mass, and virial radius of the halo, respectively, then the 
escape speed v csc is 
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2.6 km s 



4. RESULTS 

The (fWhalen et al.ll2008al ) study found several possi- 
ble final fates for halos photoevaporated by a 120 M Q 
star. If the halo is very diffuse, with n c < 1 cm~ 3 , it is 
completely evaporated anywhere it resides in the cluster, 
with no star formation. If the halo is more evolved, with 
n c > 1000 cm~ 3 , its core is shielded from both ionizing 
and LW flux and star formation proceeds there without 
delay. In halos of medium central density the relic ioniza- 
tion front shock either compresses and enriches the core 
with H2, potentially hastening its cooling and collapse 
into a star, or its residual momentum disrupts the core, 
delaying or suppressing star formation. For comparison 
to our 25 - 80 M models, whi ch are evolved to 10 Myr, 
we ran the 120 M Q models of (jWhalen et al.ll2008aft out 
to the same time. We list final outcomes for star for- 
mation in these simulations according to the criteria set 
forth in section 13.11 in Figure [2J 

4.1. 25 - 80 M G Stars 

In Figures|3]andH]we tabulate the effects of local ioniz- 
ing and LW radiation on star formation in the halo near 
25, 40, 60 and 80 M Q stars. First, we note that in each of 
these stars new star formation obeys the same trends as 
in a 120 M Q star: it is first suppressed, then delayed, and 
finally unaffected as its central baryon density and dis- 
tance to the star grow. These trends are primarily due to 



Fig. 2. — Star formation in a 1.35 X 10 s M CT ha lo in the vicin- 
ity of a 120 M star, from IWhalen et alj J2008al 1. The central 
baryon densities n c of the halo at the time of illumination were 
1.43, 10.5, 108, and 1596 cm -1 , respectively. The 120 M Q star 
was at 150, 250, 500 and 1000 pc. Completely evaporated halos 
with no star formation are labeled by crosses and halos with de- 
layed or undisturbed star formation are represented by triangles 
and circles, respectively. 



the momentum with which the I-front shock reaches the 
core of the halo. As explained earlier, when star forma- 
tion is delayed it is due either to disruption of the center 
of the halo by the relic I-front shock, LW photodissoci- 
ation of the core of the halo, or both. If star formation 
is unaffected it is because the core is completely shielded 
from LW flux, which also guarantees that the relic I-front 
shock will not reach the core before it forms a star. From 
the upper left to lower right in the 60 and 80 M panels 
runs a line that marks the boundary between quenched 
star formation and delayed or undisturbed star forma- 
tion (corresponding lines for 25 and 40 M Q stars, which 
die in SN explosions whose kinetic feedback must also be 
taken into account, are shown in section 4.2). The sup- 
pression of new star formation in the models lying just 
below this line is usually unambiguous: the shock com- 
pletely sweeps baryons from the core of the halo. Only 
in models lying just above this line does star formation 
vary with the mass of the illuminating star. The line ad- 
vances gradually up and to the right with neighbor star 
mass, but the shift is minor from 25 - 80 M Q . In a halo 
at a fixed central baryon density and distance from the 
star, the formation of a new star is at most delayed if it 
was once unaffected, or halted if it was once delayed, as 
the mass of the illuminating star varies from 25 - 80 M . 
Star formation remains unchanged above and below this 
narrow band. We note that the lines are drawn to cut 
the space between symbols evenly and as such should be 
taken as order of magnitude estimates. 
Why is star formation in the satellite halo relatively 
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Fig. 3. — Star formation in the halo when illuminated by a 25 Mq star (left) and a 40 Mq star (right) at four central baryon densities 
n c , 1.43, 10.5, 108, and 1596 cm~ 3 , and at four distances from the star, 150, 250, 500 and 1000 pc, which are typical intcrhalo spacings 
within a cluster. Completely evaporated halos with no star formation are labeled by crosses, and halos with delayed or undisturbed star 
formation are indicated by triangles and circles, respectively. 
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Fig. 4. — Star formation in the halo in the vicinity of a 60 Mq star (left) and an 80 Mq star (right) at the central gas densities and 
distances to the star considered in our study. Completely ionized halos with no star formation are labeled by crosses, and halos with delayed 
or undisturbed star formation are indicated by triangles and circles, respectively. The dotted lines mark the threshold for star formation 
in the evaporated halos, above which it proceeds and below which it is quenched. 
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uniform with neighbor star mass? On one hand, the 
dimmer flux of low-mass Pop III stars causes the front 
to transform from R-type to D-type further out from 
the center of the halo. We plot velocity profiles for the I- 
fronts at the moment they become D-type in the n c = 108 
cm~ 3 halo 500 pc from 25 - 120 M Q stars in panel (a) of 
Figure [5] In each profile the transition distance from the 
core is marked by the position of the forward peak, which 
is at 80, 75, 60, 55 and 50 pc for the 25, 40, 60, 80 and 120 
M Q stars, respectively. On the other hand, lower-mass 
stars are also longer lived, so they drive the I-front into 
the halo for greater times. This collapses the dispersion 
in I-front positions in the halo to 10 pc by the time the 
star dies, as shown in panel (b) of Figure [3 In this halo, 
the I-fronts reach the core with nearly the same velocity, 
which is well below the escape speed, with each star. 
This pattern holds for all halos in which the I-front falls 
short of the core when the star dies, when star formation 
is delayed or uninterrupted. Note that at early times the 
velocity profile is split into two smaller peaks. This is 
due to penetration of hard UV photons into the dense 
shocked gas ahead of the front, which drives a backflow 
in the frame of the shock. This is a common feature of 
I-fronts due to hard UV spectr a, as discussed in detail 
in section 4.1 and Figure 17 of llliev et all (|2009ft . The 
evolution in spectral profile from 25 - 120 M Q , which 
causes ionized gas temperatures to rise by more than 
50% in the H II region, accounts for the variation in 
peak velocity in Figure [5] Had each of these I-fronts been 
driven by a monochromatic flux of the same magnitude 
and duration, the spread in the peaks in radius would 
have bee n even less than 10 pc . 

In the iWhalen et all ([2008a) survey, the momentum 
imparted by the relic I-front shock to the core of the halo 
is primarily what determines if a new star forms there. 
Here, in a few cases, the compresssion of the shadow of 
the halo toward the axis also determines if a star forms. 
In our new models the halo forms a shadow on the same 
time scale as a 120 M star but pressure from the sur- 
rounding relic H II region drives it inward toward the axis 
for up to twice the time before the star dies, squeezing 
a flow backward into the center of the halo. This back- 
flow can slightly displace baryons from the core, as we 
show in Figure [S] for halo 073 150 pc from a 60 M Q star. 
The clump of gas originally centered in the dark matter 
potential retains its shape and is at well below the es- 
cape speed of the halo. However, its slight dislocation 
from the center of the dark matter potential at 10 Myr 
prevents it from collapsing into the star that would have 
formed if only the relic I-front shock had been present. 
Backflows are especially collimated in the axial symme- 
try we assume in our models but we still expect them 
to be present in three dimensional halos because of their 
approximate sphericity. 

In halos where stars form without interruption dur- 
ing photoevaporation, molecular hydrogen mass fractions 
rapidly rise from the initialized value of 2 xlO -6 to 10 -4 
when the star is turned on, even with LW flux, because 
the core self-shields, as we show for the halo at n c =1596 
cm" 3 500 pc from 25, 40, 60, and 80 M stars in the 
right panel of Figure [7J The core collapses even as the 
outer layers of the halo are stripped away by supersonic 
outflows, and a star forms just as quickly as in the ab- 
sence of radiation. Since we begin with cosmic mean H2 



fractions of 2 xlO -6 instead of more realistic values of 
10~ 4 for simplicity, ours are lower limits to self-shielding 
and cooling. 

One departure of the 120 M Q models from the others 
is the delayed star formation at 150 and 250 pc in the 
halo at n c = 108 cm~ 3 , which does not occur near 25 - 
80 M Q stars. This happens because the shock remnant 
traverses the core of the halo sooner, allowing gas to later 
pool in the dark matter potential and reach its original 
density at the center by 10 Myr, as shown by the dashed 
line in the left panel of Figure [7] at 150 pc. Such backfill 
allows a new star to form before a total of 20 Myr has 
elapsed, or about a merger time at z ~ 20. Similar flows 
would occur in low-mass runs evolved beyond 10 Myr but 
would probably not result in a new star prior to severe 
disruption by a merger. The higher densities at earlier 
times in the outer regions of the halo (solid lines in Figure 
[7]) are due to compression by the relic I-front shock as it 
envelopes the halo. 

We find that the higher ratio of LW to ionizing pho- 
tons of low-mass Population III stars has no effect on 
nearby star formation. The lower fluxes of these stars 
offset their higher ratios, and in any event LW photons 
from a single star cannot prevent a new star from form- 
ing in a nearby halo prior to the appearance of LW back- 
grounds at lower redshifts. They can only delay it for 
the life of the proximate star. Without LW photons halo 
collapse times at n r . ~ 2000 cm" 3 a re 7-10 Myr. At 
such densities, IWhalen et al.l ()2008al ) find that the core 
begins to strongly self-shield from local LW flux. If the 
core could collapse before the arrival of the I-front, its n c 
would be far greater, be fully shielded from LW flux, and 
form a star anyway. Halos that are completely ionized 
are too diffuse to form a star before the I-front reaches 
the core even in the absence of LW photons. Cores that 
form a star after the death of a nearby star could not 
have created one sooner without LW flux because their 
star formation times are even greater than 7 Myr. Fi- 
nally, halos that are too severely disrupted by the relic 
I-front shock to form a star cannot collapse prior to the 
arrival of the shock in the absence of LW radiation for 
the same reason. Thus, in none of these outcomes would 
a star have formed if there had been no LW flux. In the 
latter two cases, we find that H2 fractions at the center 
of the halo reacquire their original values 100 - 200 kyr 
after the star dies and that cooling and collapse of the 
core begins anew, augmented in some cases by molecu- 
lar hydrogen advected into it by the relic I-front shock. 
Any suppression or delay of star formation is entirely 
due to bulk flows driven into the core driven by relic I- 
front shock and shadow dynamics, not to the destruction 
of H2, which quickly reconstitutes in the core after the 
star dies. Thus, local ionizing UV flux governs new star 
formation in clustered halos, not local LW photons. 

An exception to this is when the baryons are very 
close to the star, 25 - 150 pc. In such cases, when 
a sta r is irradiating a clum p of baryons in the same 
halo, iHascgaw a et al.l (|2009t ) find that clouds with free- 
fall times that are shorter than ionization times can be 
prevented from collapsing by 25 M Q stars because they 
have the highest LW/ionizing photon ratio. This never 
occurs in our models because of the much lower LW fluxes 
at typical halo distances within the cluster, but such sce- 
narios are quite relevant to the formation of Population 




Fig. 5. — I-front gas velocity profiles through the central axis of the halo for 25, 40, 60, 80, and 120 Mg stars 500 pc from halo 059 (n c = 
108 cm" 3 ). Panel (a): velocity at the time each front transforms from R-type to D-type. Panel (b): gas velocity profiles at the time each 
star dies. 




Fig. 6.— Evaporated halo with n c = 1596 cm -3 (halo 073) 150 
pc from a 60 Mg star at 10 Myr. The core of the halo is slightly dis- 
placed to the left of center by backflow from the collapsed shadow 
on the right. 



Ill binaries within a halo (jTurk et al.l 120091 ). Persistent 
LW backgrounds at lower redshifts can by themselves 
prevent secondary star formation in the cluster by not al- 
lowing partially evaporated cores to cool after the death 
of the star. However, the exclusion of radiation hydro- 
dynamical effects in past studies may have led them to 
overestimate this effect, as we discuss in 5 5. 



4.2. Kinetic Feedback by SNe 

120 M stars assumed in most local 
50 M^stars die in SN ex- 



Unlike the 100 
radiative feedback studies, 25 

plosions (|Heeer fc Wooslevll2002t iTominaga et~aT]|2007[) 
Their remnants can ram-pressure strip baryons from 
nearby halos and prevent them from forming a star when 



one ot herwise would have been created. iSakuma fc Susal 
(2009) examined kinetic feedback by SNe on star for- 
mation in satellite halos with semi-analytical arguments 
that neglected radiative preprocessing of the halo by the 
progenitor star. They found that SNe generally preempt 
any star formation that fails to occur prior to the a rrival 
of the remnant at the halo. iCen fc Riquelmel (|2008j ) sim- 
ulated the interaction of SN remnants with much more 
massive halos at lower redshifts to assess the degree to 
which metals become mixed with gas deep in their inte- 
riors. These calculations, which also excluded prior pho- 
toevaporation by UV backgrounds, found some mixing 
in the outer layers of the halos due to Kelvin-Helmholtz 
instabilities but that no metals reached their interiors. 
These models did not address SNe feedback on star for- 
mation because the potential wells of the halos were much 
deeper than those in which the first stars formed, so the 
remnants could not strip gas from their cores. Detailed 
numerical models of the collision of primordial supernova 
remnants with the relic H II regions enveloping partially 
exposed cores are needed to establish the actual fate of 
star formation in these circumstances. 

Supernovae are not invariably fatal to local star for- 
mation because the time scales on which the shock 
reaches outlying halos may be longer than those on which 
baryons collapse within th em, and in s o me ca ses the rem- 
nant never reaches them. iGreif et"ail (|2007f ) found that 
SN remnants typically expand to only half the radius of 
the relic H II region of the progenitor because they come 
into pressure equil ibrium relative l y quick ly in the warm 
recombining gas. IWhalen et al.l ()200 8b) computed the 
radius of the remnant as a function of time for core col- 
lapse SNe, hypernovae, and pair-instability SNe (PISN) 
in the fossil H II regions of 10 5 - 10 7 M Q . In their simula- 
tions, ejecta from the explosions of 15 - 40 M Q stars prop- 
agate at most 400 - 500 pc in 10 Myr. We plot the radii 
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Fig. 7. — Left panel: the flow of baryons back into the center of the dark matter potential from the relic H II region, as shown in these 
density profiles along the central axis of the halo. Solid: 5.5 Myr; dashed: 10 Myr. Right panel: H2 mass fractions at the center of halo 
073 500 pc from 25, 40, 60 and 80 M stars from zero to 10 Myr. 

if the collapse time of the halo at the time of the explo- 
sion is known, Figure [3] can be modified to include both 
radiative and kinetic feedback by the star. This is done 
simply by assuming that if a star can form in the evapo- 
rated halo in the absence of an explosion, it will still form 
if central baryon collapse takes less time than required for 
the shock to reach the halo. The contraction of the cloud 
from n c ~ 10 - 10 3 t o 10 8 cm~ 3 is leisurely, taking place 
over 5 - 10 Myr (e.g. lAbel etaHl2002HO'Snea fc Normanl 
|2007| ). Thereafter, as three-body H2 cooling takes over 
in the core of the halo, the formation of a star proceeds 
very rapidly, in under a Myr. We therefore take 8 Myr as 
the average time for a partially evaporated halo to form 
a star. Clearly, this estimate is approximate because it 
neglects the initial LW dissociation of the core, but it is 
reasonable given that our studies show that H2 reforms 
there 100 - 200 kyr after the death of the star. 

Arrival times for the shock at the center of the halo 
are complicated by the fact that they are a function of 
both the energy of the explosion and the mass of its host 
halo, as shown in Figure [8l We tabulate minimum and 
maximum arrival times for the remnant at 150, 250 and 
500 pc in Tab le [D using the 15 My SN and 40 M hy- 
pernova from (jWhalen et al.l l2008b) as proxies for the 25 
and 40 Mq SNe in our study. The evolution of the 15 
Mq remnant is a reasonable approximation to that of 
the 25 Mq remnant, given the uncertainty in explosion 
energy in primordial stars over this mass range. The hy- 
pernova, whose explosive yield is ten times greater than 
that of a 40 Mq core-collapse SN, gives a lower limit to 
the arrival times of the remnant to the halo. The dis- 
persion in times with host halo mass is due to the extra 
baryons the remnant must sweep aside to reach the halo, 
which can vary by more than a factor of ten in mass. 
The host halos we have chosen, 6.9 x 10 5 M Q and 2.1 x 
10 6 Mq, bracket those in which the remnant actually es- 




10 3 
time (yr) 



Fig. 8. — Radii of 15 Mq supernova and 40 Mq hypernova rem- 
nants in 6.9 X 10 5 Mq and 2.1 X 10 6 M q halos (halol and halo2, 
respectively) from Whalcn et al. (2008b). 



of 15 and 40 M Q SN remnan ts in those models in F igure 
El As shown in Figure 10 of iWhalen et all (|2008bf ). the 
growth of these remnants does not reproduce the ide- 
alized broken power-law radii of canonical SNe, chiefly 
because they do not exhibit self-similar behavior. Mul- 
tiple reverse shocks reverberate throughout the remnant 
over its evolution, and it later violently collides with the 
dense H II region shell formed by its progenitor (see also 
iBromm et al.M2003l: iKitavama fc Yoshidall2005D . 

If we adopt the conservative stance that a SN remnant 
halts star formation in a halo as soon as it reaches it, and 
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TABLE 1 
SN Remnant Propagation Times 



SN 


halo (M Q ) 


*l50 P c(yr) 


*250pc(yr) 


*500pc(yr) 


15 M 
15 M 
40 M 

40 Mq 


6.9E+05 
2.1E+06 
6.9E+05 
2.1E+06 


1.3E+06 
3.5E+06 
3.1E+05 
8.5E+05 


5.5E+06 
> l.OE+07 
1.1E+06 
2.4E+06 


> l.OE+07 

> l.OE+07 
8.0E+06 

> l.OE+07 



capes the halo since n either progenitor can ionize halos 
greater than 10 7 M Q (jWhalen et al.ll2008bl ). Since the 
hypernova remnant, which is the fastest, travels at most 
500 pc in 10 Myr, star formation at 500 and 1000 pc in 
Figure [3] is unchanged. 

We summarize both radiative and kinetic feedback by 
25 and 40 Mq stars in Figure UU Ram pressure stripping 
at most preempts delayed star formation in the densest 
halo ionized by the 25 Mq star at 150 and 250 pc; the 
less destructive of the 25 M Q remnants only halts star 
formation at 150 pc. On the other hand, both 40 M Q 
remnants shut down delayed star formation at n c = 1596 
cm~ 3 at 150 and 250 pc. We find that kinetic feedback 
only modifies radiative feedback in halos at 400 pc or less 
from low-mass stars. The picture is more complicated 
with much more energetic PISN, which can reach any 
halo within a typi cal cluster (r < 1000 pc) in under 10 
Myr (Figure 10 of IWhalen etaLH2008bl) . Although the 
remnant overruns the halos sooner, radiative feedback 
can also speed up baryon collapse beforehand. Whether 
or not stars form in such cases remains unclear, and is 
beyond the scope of this study. 

We note that our estimates of kinetic feedback are con- 
servative for several reasons. First, as stated earlier, we 
apply arrival times for a 40 Mq hypernova to both 25 and 
40 Mq core-collapse SNe, which have lower energies and 
velocities, and hence a shorter reach within the cluster. 
Second, we assume that the halos in our study are not 
connected by filaments of dark matter and gas, which is 
usually not the case. If both the UV and the SN shock 
must propagate along a cosmological filament, with over- 
densities of 50 or more above the cosmic mean, to reach 
the core of a nearby halo they may have considerably less 
impact on subsequent star formation there. Finally, we 
also neglect the collision of the remnant with supersonic 
backflow from the evaporated halo, which will dampen 
its impact with the core of the halo. Also, 40 - 50 M Q 
primordial stars may die in more exotic explosions such 
as hypernovae o r collapsars that are e ither asymmetri c 
or beamed (e.g. iTominaga et alJl2007t iTominagal 12009( 1 . 
However, the asphericity of such events lowers the likeli- 
hood that any one halo in a cluster would be struck by 
ejecta, so for simplicity we do not consider them here. 

4.3. Analytical Fits to Local Star Formation at High 
Redshift 

In our previous study we made the distinction between 
positive, negative, and neutral feedback within in a clus- 
ter of halos. Negative feedback meant that star formation 
is either delayed or suppressed, positive feedback meant 
that a star forms more quickly in the halo than when 
there is no radiation, and neutral feedback meant that a 
star forms as fast as it would without UV. However, what 
is ultimately of importance to early structure formation 



TABLE 2 
Analytical Fits to Radiative and 
Kinetic Feedback 



M* 


ai 


6i 


ci 


25 Mq 


3.3010 


-0.33474 


-0.14613 


40 Mq 


3.1761 


-0.49480 


0.30103 


60 Mq 


3.3010 


-0.33011 


9.6910E-02 


80 Mq 


3.3010 


-0.33084 


-0.69897 



and the rise of the first stellar populations is whether or 
not a star forms in the halo, not its exact timing, which 
varies by less than a few Myr in the scenarios we have 
investigated. Since the final result of neutral, positive, 
or no feedback is the formation of a star, the distinction 
between these outcomes is not relevant. Hence, the rules 
we formulate delineate the basic threshold for star for- 
mation in a halo exposed to UV flux and SN flows from 
nearby primordial stars. 

In Figures^] and Owe plot lines that denote the thresh- 
old radius r t h from the star above which the halo can 
form a star if it is at a central gas density n c and below 
which it cannot. The threshold has the simple form 

logwnh = ai + bi{log w n c + c x ), (10) 

where rth is in pc and n c is in cm -3 . We tabulate oi, b\ 
and ci in Table [2] for 25, 40, 60 and 80 M stars, in- 
cluding both radiative and kinetic feedback for 25 and 
40 Mq stars. As shown in Figure [21 radiative feedback 
by 120 Mq stars is even simpler: if the satellite halo has 
central densities greater than 100 cm -3 a new star will 
form anywhere in the cluster. If not, its formation is 
suppressed. 

4.4. Ionization Front Instabilities 

Longer illumination times and lower fluxes promote the 
onset of dynamical instabilities in the D-type ionization 
front as it engulfs the satellite halo, as we show for the 
073 halo 150 pc away from a 25 Mq primordial star in the 
left panel of Figure [TO] They arise because the I- front as- 
sumes a cometary shape and because the high energy tail 
of the spectrum forms H2 between the front and the dense 
shell that radiatively cools the shell (e.g. iRicotti et al.l 
20011 IWhalen fc Norman1|2008a|l . Using rigorous pertur- 
bation analysis. IWilliamsl (|2002( 1 discovered that D-type 
fronts driven by photons that are oblique to the front 
are always unstable, and that the growth rates of the 
modes rise with the angle of incidence of the photons. 
This, together with cooling of the shocked shell by H2, 
is the origin of the instabilities in the I-front enveloping 
the halo in Figure [TUT 

At early times, we find that the modes with the great- 
est amplitudes are indeed those furthest out along the arc 
of the I-front, where photons are incident to the front at 
the greatest angles. At intermediate and later times the 
perturbations grow nonlinearly and degenerate into tur- 
bulent fluid motion along the outer segments of the arc. 
The amplitudes of the modes closest to the axis of the 
halo are small and instabilities never puncture its core. 
Much more prominent perturbations have been found 
in planar I-fronts approaching spherical molecular cloud 
cores in numeric al models with effic ient radiative cool- 
ing by molecules (|Mizuta et alj |2006) . These phenomena 
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Fig. 9. — Radiative and kinetic feedback on star formation near a 25 Mq star (left) and a 40 Mq star (right). Completely evaporated 
halos with no star formation are labeled by crosses, while halos with delayed or neutral star formation are marked by triangles and circles, 
respectively. The triangle overlaid on the cross in the 25 Mq panel signifies that the halo can form a delayed star if the SN goes off in 
the 6.9 X 10 5 Mq halo but not if it is in the 2.1 X 10 6 Mq halo. The dotted lines again define the boundary for star formation in the 
evaporated halos, above which it proceeds and below which it is quenched. 




50 
■g H2 mass fraction 



Fig. 10. — Temperatures at 600 kyr (left) and H 2 mass fractions at 10 Myr (right) in halo 073 (n c = 1596 cm -3 ) 150 pc from a 25 M 
star. 



have been proposed for the origin of the "Pillars of Cre- 
ation" in the Orion Nebula, but they are different from 
those in our simulat ions. They beg in as Vishniac thin- 
shell overstabilities ([Vishniac] [1983f i caused by efficient 
molecular cooling in plane-parallel I-fronts, not curved 
ones, and later erupt int o violent instabilitie s driven by 
ionizing radiation. In the lMizuta et al.l (|2006f ) models the 
unstable modes do propagate into the molecular cloud 
core. This never occurs in our simulations because H2 
cooling is too inefficient to incite Vishniac modes. 
We point out these features because they are promi- 



nent in many of our models, but they do not affect star 
formation in the halo because they never approach its 
inner regions. Mostly, they just roil gas along the shock, 
breaking it up into clumps that can persist for up to 10 
Myr. We find that they appear when the star is 150 or 
250 pc from the halo and are most prominent with 25 
and 40 M Q stars. Fewer instabilities appear as stellar 
mass increases; they arise in only two of the 120 M© 
models and have lower amplitudes. There are two rea- 
sons for this. First, higher mass stars have greater LW 
fluxes that lower H2 cooling in the dense shell. Second, 
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larger ionizing UV fluxes result in shorter-lived cometary 
profiles in which unstable modes can develop. The arc 
is crushed downward into the shadow of the halo more 
quickly and the instabilities have less time to develop. 

5. DISCUSSION AND CONCLUSION 

We find that 25 - 120 M Q primordial stars are relatively 
uniform in their effect on new star formation within clus- 
ters of small halos at high redshifts, before the rise of 
global LW backgrounds. The evolution in spectral pro- 
file from 25 - 120 M Q has no impact on the formation 
of stars in nearby halos, which allows its removal from 
the parameter space of local radiative feedback. The 
empirical fits we have devised mark the threshold for 
star formation in satellite halos as a function of central 
baryon density, proximity to the star, and neighbor star 
mass. Although the halo in our study is the just the 
least massive one found to form a star in previous AMR 
simulations, our results can be used as upper limits to 
feedback in more massive halos. Our results imply that 
future surveys of local feedback with more massive ha- 
los can be accomplished with fewer stars, since outcomes 
for halo photoevaporation above and below the belt in 
n c and radius in which there is variability is relatively 
uniform from 25 - 80 M©. 

Radiative and kinetic feedback between minihalos is 
key to many processes in early cosmological structure for- 
matio n, such as primordial SNe event rates (jWise fc Abell 
2005), especially thos e that account for cluster bias 
(jMacIntvre et al.ll2006f ). the rise of the first stellar pop- 
ulations, the assembly of primeval galaxies, and the evo- 
lution of metagalactic LW backgrounds. Our analytical 
fits enable feedback estimates for a representative cut of 
Population III stars in analytical models of these early 
processes. They can also be used in numerical simula- 
tions, especially those performed in large cosmological 
boxes capable of resolving minihalo clustering but not of 
capturing halo photoevaporation. 

Local radiative feedback at slightly lower redshifts is 
different due to the presence of LW backgrounds from pri- 
mordial stars, which until recently has been thought to be 
quite destructiv e to new star formation mediated by H2 
coolin g in halos (jHaiman et al.lll997l 120001 : lYoshida et al.l 
2003). However, recent, more detailed simulations re- 
veal that star formation in cosmological halos is post- 
poned rather than prevented in LW backgrounds, even 
large one s that are consistent with a fully re i onize d 
universe (IWise fc Abell [20071 : IfTShea fc Normanl l2008h . 
IWise fc Abell (120071 ) found that a halo that formed a star 
at 5 x 10 5 M© in the absence of a photodissociative back- 
ground still formed one by H2 cooling 50 Myr later after 
it grew by mergers and accretion to 5 x 10 6 M© in a 
uniform LW field of 1 J 2 i (= 10~ 21 erg cm~ 2 Hz -1 str -1 
s _1 ), that of a fully-ionized universe. Other halos in its 
vicinity also grew to larger masses, even though cool- 
ing and collapse of baryons were temporarily stalled in 
them. Thus, at lower redshifts local UV feedback still be- 
gins when one star irradiates neighbor halos. In contrast 
to the first generation, LW backgrounds may continue to 
suppress star formation in photoevaporated halos after 
the death of the star by slowing the formation of H2 in 



the relic H II region, the remnant I-front shock, and the 
halo core. 

However, this effect may have been overestimated in 
previous analyses. Consider the morphology and large 
densities and H2 fractions in the relic H II region at 10 
Myr in the right panel of Figure 1101 Similar H2 fractions 
would likely persist in the recombining H II region even 
in large LW backgrounds because the very high electron 
fractions there restore it so quickly via the H~ channel. 
This is especially true in the high densities of the I-front 
shock remnant, that can be seen to envelope the core 
in Figure [TOl Enough molecular hydrogen could be cat- 
alyzed in the envelope to expel the LW background from 
its interior, allowing H2 to reform at the center of the 
halo, cool it, and form a star. Thus, photoevaporation 
may actually free satellite halos to form stars that were 
previously suppressed by the background. This, together 
with our current study, suggests that star formation in 
cosmological halos at lower redshift was much more ro- 
bust than is often supposed. Numerical models are now 
being developed to investigate the survival of molecular 
hydrogen, and therefore new star formation, in evapo- 
rated halos for a range LW backgrounds. 

The effect of photoevaporation on the final mass of any 
star that does form in t he halo is not yet w ell understood, 
but initial estimates bv lSusa et al.l ()2009l ) suggest that it 
will be smaller than in undisturbed halos. They find that 
outflows and shock disruption in the core lower central 
accretion rates, and by extrapolating these rates from 
early stages of collapse out to Kelvin-Hclmholtz contrac- 
tion time scales, they conclude that the final star will 
be 25 - 50 M©. This mass scale is similar to those on 
which Pop I II.2 stars form due to HD cooling in relic 
H II regions (Yoshida et al. 2007). HD is important be- 
cause it can cool primordial H II regions down to the 
CMB temperature and lower the mass scales on which 
they fragment. However, we do not include it in our 
models because it forms primarily in the relic ionized 
gas surrounding the halo core, not in the core itself. 

Although our study is a significant extension to our 
earlier survey of local radiative feedback, additional feed- 
back channels remain to be properly investigated. If 
the death of the star results in a black hole, accre- 
tion would expose nearby halos to its s oft x-ray flux 
(jMachacek et al.ll2003t [Alvarez et al.ll2009f h creating sig- 
nificant free electron fractions in them due to secondary 
ionizations without strongly heating them. This process 
could enhance their H2 mass fractions and promote their 
collapse into new stars. Likewise, the impact of SN ejecta 
with a halo that has been partially stripped by supersonic 
flows could deposit metals into its interior and accelerate 
its cooling and collapse. These potential avenues of pos- 
itive feedback on primordial star formation will be the 
focus of future simulations. 



This work was carried out under the auspices of the 
National Nuclear Security Administration of the U.S. De- 
partment of Energy at Los Alamos National Laboratory 
under Contract No. DE-AC52-06NA25396. The simula- 
tions were performed on the open cluster Coyote at Los 
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